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Abstract 

Using some modification of the standard fermion technique we derive factorized formula for 
spin operator matrix elements (form-factors) between general eigenstates of the Hamiltonian 
of quantum Ising chain in a transverse field of finite length. The derivation is based on the 
approach recently used to derive factorized formula for Zj\r-spin operator matrix elements be- 
tween ground eigenstates of the Hamiltonian of the Z/v-symmetric superintegrable chiral Potts 
quantum chain. The obtained factorized formulas for the matrix elements of Ising chain coincide 
with the corresponding expressions obtained by the Separation of Variables Method. 
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1 Introduction 

During the last decade more attention has been drawn to spin matrix elements in finite Ising 
systems. Bugrij and Lisovyy in [JJ[2] has been pointed out that one may write completely factor- 
ized closed expressions for spin matrix elements and form-factor representation for correlation 
function of the two-dimensional Ising model on finite lattice. The finite volume form-factors 
of the spin field in the Ising field theory was derived by Fonseca and Zamolodchikov [3]. The 
Bugrij-Lisovyy formula for spin matrix elements on finite lattice was proved in [H [5] by the 
Separation of Variables Method. In [U [7] another approach to calculation of spin matrix ele- 
ments was suggested. While it allows to get the answer in the continuum limit on the cylinder, 
the lattice situation appears to be more problematic. 

Using the Separation of Variables Method the factorized formulas for spin operator matrix 
elements were obtained for the quantum Ising model in a transverse field [H [5] and for the 
XY quantum chain [8]. The quantum Ising chain (QIC) is the special case (N = 2) of the 
Zjv-symmetric superintegrable chiral Potts quantum chain (SCPC) [9]. In the framework of 
extended Onsager algebra [1QJ Baxter conjectured the SCPC-spin matrix elements for the ground 
state Onsager sectors |11J . Independently in the same framework in [12] these matrix elements 
were re-obtained, generalized to arbitrary Onsager sectors and proved. There it was shown that 
the matrix elements of spin operators between the eigenvectors of the SCPC Hamiltonian for 
arbitrary Onsager sectors can be presented in a factorized form with unknown normalization 
factors AfpQ depending only on Onsager sectors but not on particular vectors of these sectors. 
The derivation does not use the information on how these Onsager sectors are included in 
large quantum space of the model. But in order to calculate scalar factors A/pg we need such 
information. Thus probably the derivation of these factors will require an information on the 
Bethe-states of related T2-model and sZ(2)-loop algebra symmetries [TB I f!4" l [T5l [To] IT7]. 
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In |10| Baxter noted that the calculation of the matrix elements of spin operators between 
eigenvectors of Hamiltonian H = Hq + hH\ is greatly simplified if one starts from the calculation 
of matrix elements of spin operators between eigenvectors of Hq (that is of H at h = 0). In this 
paper we use fermion technique to describe the Onsager sectors of QIC. To identify the bases of 
Onsager sectors in [12] with fermion bases we diagonalize Hq by means of dual Jordan- Wigner 
fermion operators. As in the case of the duality transformation of 2D Ising model on finite lattice 
|18] , the duality transformation for the finite quantum Ising chain requires consideration of the 
periodic and antiperiodic boundary conditions in a uniform way together with a modification 
of the standard Jordan-Wigner fermions. Then derivation of spin matrix elements between 
eigenstates of Hq reduces to Cauchy determinant which in turn allows to extract the factors 
Npq = Np 0t Q in an explicit form. In order to get matrix elements of spin operator between 
eigenstates of H we use Bogoliubov transformation of fermion operators. This transformation 
of a pair of fermion operators with opposite momenta is directly related to rotation in a two- 
dimensional tensor component of Onsager sectors. Bogoliubov transformation leads to a sum 
over the different pairs of fermions which can be summed to get a factorized formula (180p for 
matrix elements of spin operator between eigenvectors of H. 

This paper is organized as follows: In Sect. 2 we define the Hamiltonian of QIC with unified 
periodic and antiperiodic boundary conditions, duality transformation, Jordan-Wigner and dual 
Jordan-Wigner fermions. In Sect. 3 the derivation of matrix elements of spin operator between 
eigenstates of Hq is presented by reducing them to Cauchy determinant. This allows to find the 
scalar factors Np 0i Q a (Appendix A). In Sect. 4 we define Bogoliubov transformation of fermions 
and find relation between eigenstates of Hq and H\. In Sect. 5 we use the results of Sects. 3 and 
4 to get factorized matrix elements of spin operator between eigenstates of H. The derivation 
uses summation formula from Appendix B. Sect. 6 summarizes our results. 

2 The quantum Ising chain 

2.1 The Hamiltonian of the quantum Ising chain 

The quantum Ising chain (QIC) of length L with periodic boundary condition is defined by the 
Hamiltonian 

L-1 L 

H p = fljf + hH[ = - £ a%a x k+l - a x L af - h £ a% . (1) 

k=l k=l 

The same model with antiperiodic boundary condition is governed by the Hamiltonian 

L-1 L 

H A = H$ + hH? = -J2 o%a% +1 + of of - h °l ■ (2) 

k=l k=l 

Both Hamiltonians include term describing interaction with a transverse magnetic field of 
strength h . The space of states of these systems is the L-fold tensor product of two-dimensional 
spaces Vj, j = 1, . . . , L: V = Vi<8> • • • ®Vl- The spin operators 

of = 1®...®1 ® v o^®l <S>... ® 1, 

k-th 

where a = x,y, z, k = 1, . . . , L, and a a are Pauli matrices: 
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act non-trivially only on the space Vk in the tensor product V. In order to calculate matrix 
elements of the spin operator a k in this model it is convenient to consider the periodic and 
antiperiodic boundary conditions in a unified formalism. For this aim we introduce additional 
auxiliary two-dimensional space Vb- The space of states becomes 



and unified Hamiltonian is 
where 

L-l 



V = FfeOViO • • • ®V L (3) 
H = H + hfli, (4) 



H o = -^2^k a k+i -^l^i = -^2°k a k+ii H i = -^2a z k , (5) 

k=l k=l k=l 

with the following boundary conditions for spin operators 

(T L+ i = cr b a l , a L+1 =a 1 , (6) 
which can be extended to a "quasiperiodic" conditions 

a k+L = (J b (T ki cr k+L = cr k- {') 

In what follows the main object of our consideration is the Hamiltonian Q acting on the 
space (J3]). The space of states V decomposes into direct sum V = V P ®V A of two subspaces 
V p and V A with eigenvalues +1 and —1 of cr|, respectively. Since <t| commutes with the 
Hamiltonian H these two subspaces are invariant under the action of H. The restriction of H 
to these subspaces is H p and H A , respectively. The boundary conditions © become periodic 
conditions, = af, & z L+l = af, or antiperiodic condition, <r£ +1 = —erf, (J z L+l = erf, with 

respect to the eigenvalue of a^. 

The Hamiltonian H also commutes with the Hermitian operator W 

W = -at f[ ol . (8) 

fc=i 

Since W 2 = 1, the eigenvalues of W are w = ±1 and the space H can be decomposed into 
the direct sum V = Vns®Vr of two subspaces Vns and Vr corresponding to the eigenvalues 
w = — 1 and w = +1, respectively. The subspace Vns (resp. Vr) is called Neveu-Schwarz 
sector or NS -sector (resp. Ramond sector or i?-sector). 

With respect to the eigenvalues of the commuting operators and W, the space of states 
V decomposes into the following direct sum: 

V = V A S ®V P S ®V A ®V P , (9) 

where, for example, if |^) € V NS we have Wj^) = — \^>), cf| x I / ) = — Each of these subspaces 
have dimension 2 L_1 and is invariant under the action of H, Hq and Hi because they commute 
with a% and W. 

Let us define the translation operator on V 

T = \ (1 + a% + of - o z h a\) T 1>2 T 2>3 • • • T L -i,l , (10) 



where the operator 



Tfc, fc+ i = - (1 + a x k a x k+1 + a y k a v k+1 + a z k a z k+1 ) 
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swaps the vectors in the spaces Vj~ and Vk+i- The translation operator T acts on the local spin 
operators as (see also flS}) 

T^ = ^ +1 T, Tat = at +1 J, Taf = a z b T, Tof = of erf T (11) 

and satisfies the relations 

TT+ = 1, T L = - (1 + a$ + W - o%W) , J 2L = 1, TW = WT . 

QIC possesses generalized translation invariance with respect to T, i.e. [H, T] = 0. On the 
subspace V p the operator T turns into standard translation operator. 

2.2 Duality transformation 

The presence of the matrix <r| in the Hamiltonian results in some modification of standard 
duality transformation of QIC. We place the dual Pauli matrices af, I = 1, . . ., L, on the dual 
chain sites which are located between the neighboring sites of the original chain, and take the 
following numbering of dual sites: the dual site I is located between the sites I and I + 1 of the 
original chain. We define the dual Pauli matrices by the following relations: 

l L 

~X X I I Z ~Z X X ~X X Z ~X X I I z 

a i = a b 11 °fc> a l = a l a l+l' cr l =(7 b< J l> cr i = cr 6ll cr n 

k=l 1=1 

L 

^ Z ~ Z ^r X ~ X 7^ Z ^ X Zr X * Z ~ X ~ 2 I I ~ Z /lO\ 

°L = °b a l a Li <7l=<7l<T2i cr 6= cr 6 CJ l> (J b=[[ cr l- I 12 ) 

1=1 

The inverse duality transformation has the form 

L L 
„x ~x I I ~z z ~x ~x x ~x ~z x ~x I I ~z 

a k = a b . a k = a k-l a ki (T L = cr b cr L> ®1 = ®b 11'' 

l=k 1=1 

L 

„Z ~X ~X Z ~Z~X~X X ~Z ~X Z I I ~Z /1Q\ 

a L = a L-l a L, <J l= (J b a l (J L^ °b= (J b< J Li (J b=[[ cr l- \ 1A ) 

1=1 

Duality transformation interchange Hq and H\ and in terms of the dual Pauli matrices the 
Hamiltonian H becomes 

H = H + h#i = h (- Sfdf+i ~ h Yl J = h + hHi) = h#, (14) 

V i=i fc=i / 

where h = h" 1 and H is the dual Hamiltonian with boundary conditions 

~x ~z~a; ~z ~z f-\K\ 

a L+1 = a b a 1 , a L+1 = a 1 . (15) 
Let us define the dual operator W by 

L 



i=i 

It satisfies the relation W = W, that is W is invariant under the duality transformation. 
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2.3 Jordan— Wigner fermion operators 

Usually the first step for diagonalization of QIC Hamiltonian is the introduction of the Jordan- 
Wigner (JW) fermion creation and annihilation operators. In our case they are 

fe-i fe-i 

a* = °i n °i • a k> 4 = a b n a i • > 

1=1 i=i 

where = (cr^ ± ia v k )/2. It is easy to verify that JW operators (fT6|) satisfy: 

a) the standard anticommutation relations 

{c k ,cf} = S k ,i , {cfe,c/} = 0, {c^,c / + } = 0, (17) 

b) the boundary conditions (we take into account © and cr£ , 1 = cr^af) 

CL+1 =\Nc u c + +1 = Wc+, (18) 

c) the "quasiperiodic" conditions 

c k+L = \Nc k , c + +i = Wc+, (19) 

d) action of the translation operator T consistent with (|19p 

Tcfc = Cfc+iT, Tc^ = c^ +1 T . (20) 

In terms of the fermion number operator Q = ^fe=i c k c k we have the relations 

L L 

W = -<x 6 2 [J ^ = II exp(ivrc+c fc ) = (21) 

fc=l k=l 

[W,c fc ] = 0, [W,c+]=0. (22) 
The operators -Ho an( A f l from ((5|) in terms of the JW fermion operators take the form 

L L 

Ho = -Yl ^t+i = - £(4 " c ^)(4+i + <%+!)> (23) 

fc=i fc=i 

Hi = -^ = 2 J2 U^k ~l)= E(4 - + c%). (24) 

fc=i fc=i ^ ' fc=i 

The last equation gives that JW fermions diagonalize Hi. 



2.4 Dual Jordan— Wigner fermion operators 

As discussed in the Introduction in order to calculate the spin operator matrix elements in the 
QIC we must first solve the auxiliary problem, namely, to calculate the matrix element in the 
basis of eigenvectors of the Hamiltonian Hq, which is H at h = 0. In order to diagonalize 
Hq it is convenient to define the dual JW fermion creation and annihilation operators at, a k , 
k = 1, . . . , L, through the relations 

c k ~ c k = ~( a k ~ Ofe), c£ +1 + c k+ i = a+ +a k (25) 
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interchanging Hq and H\ given by (|23j) and (|24|) : 

L L 



fc=l fe=l fc=i 

L L 

Hi = ^2(4 - c k){4 +<*) = - X]( a fc+i - + a fc) • (27) 

fe=l fc=i 

From (|25p. we have the explicit formulas for the dual JW fermion operators 

ak = ^ (4+1 + c ^+i + ( c fc ~ c ^)) ' a fc = 2 ( c fc+i + Cfe+1 ~ ^ c fc ~ Cfc )) ' ( 28 ) 
Also these operators can be written in terms of a k 

k 

ak 



\°% n <* - <*) > 4 = \°t ■ n < • K+i + • ( 2g ) 



2 

z=i «=i 



or in terms of the dual operators crj? 



L L 
Z=fc+1 Z=fc+1 



where a k = (cfj£ ± ia y k )/2 and <t| = icr%crl ■ 

Due to the linearity of the transformation (|28|) and using the properties (|17p ~ (|22|) of the 
JW operators cf, q together with (j29|) . a straightforward calculation gives for the dual JW 
operators: 

a) the anticommutation relations 

{a k , af} = 8 kj i , {a k ,ai} = 0, {a+,a+} = 0, (31) 

b) the boundary conditions 

a L+1 = Wai , a+ +1 = Wa+ , (32) 

c) the "quasiperiodic" conditions 

a k+L = \Na k , a~l +L = Wa^ , (33) 

d) the commutativity with the operator W = W 

[\N,a k }= [W,a+] =0, (34) 

e) the anticommutativity with the operator of 

K,a fc } = 0, K,a+} = 0, (35) 

f) action of the translation operator T consistent with (|33p 

Ta k = a k+1 J, Ta~l = a^ +1 J . (36) 



From (|34p it follows that the subspaces Vns an d V/j are invariant with respect to the action 
of fermion algebra ([31 P having unique irreducible representation — standard Fock representation 
of dimension 2 L . Since the dimensions of Vns an d Vr are also 2 L , the action of the algebra 
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(|3ip on each of them is irreducible. We will denote the corresponding vacuum states in the 
irreducible subspaces Vns and Vr by |0)o,at5 and |0)o,_r, respectively. Using ((5]) and (J22D it is 
easy to verify that the following two vectors 



|0) 



1 



0,NS 



2 (L+l)/2 



(6) 



(1) 



+ 



(6) 



(1) 



\0)0,R 



1 



(6) 



belonging to V 
states: 



Also we have T|0) 



2 (L+l)/2 

v N s®v R 



(i) 



(L) 



(6) 



1 



(1) 



1 



(37) 
(38) 



W|0) 



0,NS 



■|0)o, 



W|0) ,fl = |0) ,fl, 



®Vl satisfy the necessary equations for these vacuum 

ai\0)o,NS = 0, l = l,...,L, 
ai\0) o ,R = 0, 1 = 1,. ..,L. 



0,NS 



|0)o,tvs and T|0) p = |0) ,r 



Because of ()34|) we may restrict operators a;, , I = 1, . . . , L, to the Neveu-Schwarz sector 
Vats (resp. to the Ramond sector V/j) where W has eigenvalue —1 (resp. +1). The restricted 
operators will be denoted //, (resp. di, d^~) which have due to ([32]) antiperiodic (resp. 
periodic) boundary condition: 



/ 



L+l 



~4 



I ' 



4 + i 



4- 



(39) 



The set of the vectors 



(ftr(m n2 ---(m nL \o) ,Ns, (^rwr-KnoKn, n fc e{ ,i}, ( 40 ) 

constitute a basis of V = Vat^ © Vr on which the operator ([26]) is diagonal. 

As it was explained in the Sect. 2.1, the boundary condition of H and iio is fixed by the 
eigenvalues of the operator a?. Using (|30p one can relate the dual fermion number operator Q 
to of: 

L L 

of = II 5 f = (- 1 ) <i ' Q = E a ^- ( 41 ) 

1=1 k=l 
Hence the states from ([40 p with an even number of excitations (including the vacuum states 
|0}o,JV5 and |0)o,j?) belong to V p and the states with odd number of excitations belong to V A . 

Let us define the momentum representation for the dual Jordan-Wigner fermion operators. 
It is defined by means of the discrete Fourier transformation: 



1 - 



d„ 



1 L 



iq I 



-ipl 



i=i 



l 



(42) 



(43) 



where in the A^S'-sector, due to the antiperiodicity ([39p . the momentum q takes "half-integral" 
values q = 27r(fe+l/2)/L, k G Z, and in the i?-sector, due to the periodicity ([39]) . the momentum 
p takes "integral" values p = 2irk/L, k G Z. All the momenta are defined up to a multiple of 
27r. The set of "half-integral" (resp. "integral") momenta q with — ir < q < tt will be denoted Q 
(resp. V). All the momenta enter Q (resp. V) by pairs (i.e. two different momenta q, —q £ Q) 
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except the momenta and ir. The momentum G P, the momentum ir G V if the length L of 
the chain is even and ir G Q if L is odd. In what follows we will use the basis 

\Q)o,NS = kl, • • • , <?m)o,7V,S = f qi f q2 " " " /^|0}o,JVS £ V7VS , (44) 

|-P)o,.R = \Pi, ■ ■ ■ ) Pn)o,R = ■ ■ • d+ n \0) OiR e Vr (45) 

of V instead of basis (j40j) . These states are labeled by the subsets Q C Q and P C P of excited 
fermion momenta in NS-sector and R-sector, respectively. We introduce the subsets Q+, Q_, 
Qo C Q as 

Q+ = {q G Q |0 < q < vr , -</ G Q}, Q_ = G Q \ - vr < q < , -q G Q}, Q = Q\(Q+ U Q-) 

and the subsets Q + , Qo C Q = Q\Q by similar formulas with the replacement Q — > Q. 
The cardinality of the sets Q, Q + , etc. are to, to + , etc., respectively. We also will use the sets 
Q + = Q + U Q+, Qo = Qo U Ooi Q- = Q- U Analogously we introduce the subsets P, P + , 
etc. of P which have the cardinality n, n + , etc. and the sets V+ = P+ U P + , Po = P) U Po, 
P_ = P_ U P_ . We have the relations 

L = m + fh = n + n, m = 2m + + too, to = 2to+ + mo, n = 2n + + no, n = 2n + + no- 

Formulas (|42p and (|43p imply that the translation operator T acts diagonally on fermion 
operators in momentum representation: 

T/ g = e J %T, T/+ = e-*V+T, Td p = e*%T, Td+ = e^d+T . (46) 

Using (|26p and (|27p we can rewrite restrictions of Po an d Pi on NS and R-sectors in terms of 
dual Jordan-Wigner fermion operators in momentum representation: 

P = ^(2a+a g - 1) , Hi = ((1 - 2a+a q ) cosq + i (a+atg + a q a_ q ) sing) , (47) 

where for iVS'-sector a q = f q , q G Q, and for P-sector a q = d q , q G P. In what follows 
we will often use the same notations H, Hq and Hi for the restrictions of the corresponding 
Hamiltonians (|4" ]) ; ([5 ]) to subspaces Vns and Vr. 

3 Spin operator matrix elements 

between the eigenstates of the Hamiltonian Hq 

In this section we will find the matrix elements of spin operator af between eigenstates (|44p 
and (I45p of Po- The relation \Naf = —af\N gives that the action of af swaps two eigenvalues 
±1 of W, that is the operator af maps Vns to Vr and vice versa. The operator af commutes 
with Po from ([5]) and therefore due to (|26p it commutes with the dual fermion number operator 
Q from (141 p . Hence action of af conserve the number of excitations. In particular, it should 
map |0)o,at5 to a vector proportional to |0)o,_r. Using explicit formulas ([37]) and (f38|) we get 

|0)o,Ji = <7?|0)o,j\ra, |0)o,jv5 = ^i|0) 0)J r. (48) 
Using (I29p we get the commutation relations 

a t af + afai = Q, afaf + afaf = 0, 1 = 1,2,... ,L, 

afd k af = -f k , afdfa x i=-f+ l = l,2,...,L, (49) 



S 



which give after Fourier transformations (|42p and ([4~3 



X 1 _x 

a 1 a p a 1 



(g-p) k 



Vl 



k q 



+ p-i{1-p) k 



q 

Z^ J i 1 _ e - i 



(50) 



(51) 



allowing to calculate matrix elements o,r(P\&i \Q)o,ns of the spin operator af between vectors 
(|44p and (|45p . Since af conserve the number of excitation, the non-zero matrix elements will 
appear only if n = m: 

o,R{Ph Pn\of\qi, q m )o,NS = oM°\ d Pn ■ ' • d P2 d pi a i /+/+ • • • f+ m \0) ,ns = 

„/(\\A A ...A rr X f r+ f+ . . . f + 

_„ 1 - e~ l 
f-1) 



2 

~ Z 



-i(«-J>7) 0,-R( Mpn d Pm-l • • • dp 2 affqf qi f q2 ■ ■ ■ /,J0) ,JVS 



J n,m 



fc=i 



-ife-pi) 



o,fl(0|dp m 4 



Pm-l 



^i7+4 • • • /S-!^ • • • £J0W (52) 



Taking into account o,-R(0|of |0)o,jvs = 1 which follows from ([IS]), we see that we get recurrent 
relation for the expansion of the determinant 



(53) 







A,2, • 




\ 


0,R(°\ d Pm ■ ■ • d P2 d Pi a lfqJq 2 ' ' ' fq\J°)o,NS = S n ,m det 


^2,1, 


4j,2, • 


^-2,m 






\ An,!' 


An,2, • 




/ 



with the matrix elements 



.4 



1 



L J e iq i - 



(54) 



with respect to first row. It is a variant of Wick theorem with two-particle pairing given by 
Aij. The determinant (|53[) can be transformed to the Cauchy determinant for which we have 



det 



1 



III 



<i<j< 



n ( x j ~ x i) (Vi ~ Vj) 



, x i Vj J X<i,j<n Yli,j=l( x i Vj) 

Introducing the variables Xj = e iqj , yi = e tpi we get 

o,r(-P|o"i \Q)q,ns = o,h(Pi, P2, ■ ■ ■ ,PmWf\qi, q2, ■ ■ ■ , q m )o,NS 



< m .(-l) m(m - 1)/2 ( ^ 



L 



m piqt 



n?A eiqi - eiPj ) 



(55) 



For convenience, we will use the following ordering of the momenta in the state |Q)o,jvs = 
\qi, q2, ■ ■ ■ , q m )o,NS : from the beginning the pairs (—q,q), for q G Q + , are going on and then 
q € Qo are going on in an order fixed over all this paper. The final formula for the matrix 
elements ([80]) will not depend on this ordering. Such ordering defines the sign of \Q)o,NS 
uniquely. Analogously we define the ordering of the momenta for \P)o,R- It is obvious that the 
sets Q + , Qo and P + , Po completely characterize the states \Q)o,ns and \P)o i r- Hence we can 
use the following notations: 



\Q+, Qo)o,NS = \Q)o,NS i \P+, Po)o,R = \P)o,R 



(56) 
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The operators Hq and H\ generate Onsager algebra. From (|47p it follows that the space Vq C 
Vns of all the states with the same fixed Qq is invariant with respect to the action of Hq and H\ 
and in fact it is space of an irreducible representation of this Onsager algebra. The basis vectors 
in Vq are labeled by the sets Q+ and the dimension of Vq is 2l 2 +L The space Vq can be 
identified with the tensor product of |Q+| two-dimensional spaces with each tensor component 
labeled by q G Q + . To vector i~~ m + \Q+, Qo)o,7VS we put in correspondence the tensor product 
of two-dimensional vectors in each tensor component: (1,0) T if q G Q+ and (0, 1) T if q G Q+. 
Then from (147p it follows that in the basis of vectors i~ m+ \Q+, Qo)o,ns with above mentioned 
identification of bases the Hamiltonian H has the form of Eq. (10) from [12]. We have analogous 
situation for Onsager sector Vp Q C Vr. 

The results of [12] predict the matrix elements of spin operator up to unknown constants 
Np 0j q q depending only on Onsager sectors: 



0,R{ P \°l\Q)0,NS = 0l R(P+> P oWl\Q+,Qo)o,NS = i {na - ma)l2 S 2 \I\,L-m -n N Po ,Q ^ 

x (~l) n+ n ae /tanf 

n ae/ (cos a + 1)° r[/3e/( cos f 3 ~ X ) r Uaei U/3eT( cos a ~ cos I 3 ) ' 



(57) 



where 



o,o- 



I = Q + UP+, I = P + UQ + , \I\=m + + n + , 
L = 2m + + 2fh + + itiq + mo, 2\I\ — L + no + mo = m — n, 
2r = -P 0>0 + P ,o + Qo,o ~ Qo,o + 1, 2ff = -Q 0)7T + Q ,^ + Po.tt - Po.tt + 1, (58) 
1, OgPq o fl, ttGPq s f 1,0GP o fl, ttGPo 



0, 0G*P ' Po H o! ^P ' Po H o! O0P O ' P °H 0, ^P ■ (59) 



The functions Qo,o, Qo,tt, Qo,o> Qo,tt are defined similarly. Note that always we have Qo,o = 
Oo,o = but we included them for uniformity. From (|58p it follows a, r G {0, 1}. Also we used 

Om,n = ^mo-noi2(n+-m+) = "2\I\,L— mo— no ' 

Formula ([55]) fixes Np 0j q in ([57]) uniquely. The derivation in Appendix A gives independent 
proof of the formula (j57[) (it is the identity f)81 1) at n = m) and allows to fix Np 0j q in different 
equivalent ways. One possible expression for Np ^Q is given by 



r(mod2) [2\ 



(mo+no)/2 



where 



NP ,Q = Co_^0 ) ^- ^2 iVe A Po)Qo X 

x n ^7 /2 ) n (cosg-i) T n ( C0S p+ir n n (^s P -co Sq ), m 

jyO _ ^_^-jno(no-l)/2^_^-j(mo-™o)(mo-n -2)/8^_^^(L--r?io-no)(no-m.o)/4 
_/y 2 — 2-( n 0-m ) 2 /4^ A" e = JJ e ip(m -n )/2 JJ e -iq(m -n -2)/2^ 

pefo geQo 

A _ U.p<p'GP ( eW ~ e%P ) Wq<q'&Q ( e%q ~ e%<1 ) r(mod2) _ 1 t , x no - mo , -, \ 

It can be obtained from (|8ip by fixing P + and Q+ as empty sets. This expression is explicitly 
does not depend on particular vectors of Onsager sectors. This expression for Np Q Q will not 
be used for the calculations in the paper. Another more useful equivalent expression for Np 0i Q 
is given by ([75]) . Now we can substitute the obtained Np 0j Q to the formulas of Sect. 4.2 and 
Sect. 4.3 of [12] (they correspond to formula (|77p of present paper) and obtain after some 
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regrouping of factors the factorized matrix elements of spin operator between eigenvectors of H. 
Instead of this direct way we choose more long but self-contained presentation of results and 
show also that the rotations in representations of Onsager algebra correspond to Bogoliubov 
transformations of pairs of fermions with opposite momenta. 



4 Eigenvectors of the quantum Ising chain 

Using (|47p we can rewrite the initial Hamiltonian H = Hq + \\H\ in terms of dual Jordan- Wigner 
fermion operators in the momentum representation: 

H = ^2 {(^ a t a g — — h cos q) + i (a~^a~^ q + a q a_ q ) h sin q) , (61) 

where for iVS-sector a q = f q , q € Q, and for i?-sector a q = d q , q £ V . Bogoliubov transformation 
defines new fermion operators, < q < tt: 

b q — C q CL q -\- is q Cl__q, bq — C q CL q is q Q>— q , 
b~ q — CqCL — q iSqQig , b_ „ — C q CL_ q + is q CL q: 



(62) 



which diagonalize the Hamiltonian (I6ip 

where in the ferromagnetic regime, < h < 1, we have 



e q = yl + h 2 — 2hcosg , £q = 1 — h , e 7r = l + h, (63) 



h sing 
1 — h cos q 

Note that c q and s q can be represented in the form 



— _£ 

Cq — COS ^ > Sq 



■ "Q 

sin — 
2 



tanfl g = - (64) 



(e w + e q )(e + e q 



c q {e q , h) = y ^ ' ^ 2_, sg(e g , h) = — - e q )(e q - e )e q . 

Formally we have 

c<z(-£ g ,h) = s g (e g ,h), s 9 (-e 9 , h) = -c g (e 9 , h) . (65) 
The vacuum states |0)jv5 an d |0)/j for the Hamiltonian in iVS 1 - and .R-sector 

\0)ns= II ( c <? + is if-qfq) \°)o,ns , \0)r= J] (c g + ^l 9 4)|0) , R , (66) 

0<g<7r 0<g<7r 

where |0)o jvs an d |0)o,i? are the vacuum states of the Hamiltonian Hq in corresponding sectors. 
It is easy to show that for both sectors the following relations are fulfilled: b q \0) =0 and 

bq {c q + is q at q a+) |0) = a+|0} 0) b± q b+ (c q + is q a± q a+) |0) = (c 9 al 9 o+ + is q )\0) . (67) 

As in the case of the eigenvectors of Hq we will label the eigenvectors of H by the same type 
of sets of momenta and use the same ordering. For example, in the NS-sector the vectors are 
labeled by a set Q of momenta: 

\Q)ns = \qi, q m )Ns = b qi b+ 2 ■ ■ ■ b+ m \0) NS g V NS 
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or by two sets and Qq with the following ordering of the momenta in the state: 



\qi, Q2,---, q m )NS = \Q)ns = \Q+,Qo)ns = {bt q b+) b+ | 0) NS ■ 

<?£<?+ q&Qo 

The formulas (|46p for the action of translation operator together with the formulas for Bogoli- 
ubov transformation imply 

Jb q = e% q J, Jb+ = e-%+T 

and therefore 

T\Q)ns = e- l ^Q q \Q) NS , ns{Q\T = Ns{Q\e- lT ">^ q . (68) 



The formulas (|66p and (|67p allow to find explicit expressions for the eigenvectors of H as 
linear combinations of the eigenvectors of Hq: 

\Q+,Qo)nS= n {Cqa- q aq + iSq) Y[ {°q + is q a -q a t) IJ a gl°)o,JVS- (69) 

<je<2+ q&Q+ q^Qo 

It is convenient to introduce the notations: e(q) = e q if q 6 Q+ and e(g) = — e q if g G Q+, 



^q) = J^^^ ± ^ 1 , /3(g) = ^yV(^-^))(^)-^oMg). (70) 

The formulas ([65]) give 

a (?) = Cg, /5(g) = s g if q£Q+; a(q) = s q , /3(g) = -c 9 if q€Q + , 
which allow to present the products over Q+ and in (169p uniformly: 



NS 



\Q+,Qo)NS=i m+ ' {a(q)a± q a+ + i(3(q)) ] [ a+|0) . 



= i m +{ _ 1) m + C™'+ J] /3(g) J] «(<?) lQ+,Qo>o,*s, (71) 

where 

|Q+,Qo)o,.ra = J| ( a -<? a ?) II a g"l )o,JV5- 
<jeQ+ <?eQo 

It means that the eigenvectors of the Hamiltonian H in given Onsager sector Vq can be 
represented as linear combinations of the eigenvectors of the Hamiltonian Hq from the same 
Onsager sector. 

Similarly to (17ip the eigenvectors of the Hamiltonian H from i?-sector are 



\P + ,P ) R = i n n-l) n + Yl rn ' + II M II a(p)l^o>o,fl, (72) 

p;cp+ pgp; pep; 

where 

m,Po)o,p= n n <i°>o,fl- 
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5 Spin operator matrix elements for the eigenstates 
of Hamiltonian H 

For calculation of the matrix element R(P\crf \Q)ns we use the expressions (|7T|) and ([72]) for the 
eigenstates of the Hamiltonian H and (|57|) for the matrix element o,r(Po, P+|crf \Q' + , Qo)o,ns- 
In order to obtain the factorized formula for the matrix elements n{P\af\Q) ms we need to make 
summation over the eigenstates of Hq labeled by the sets Q' + and Pi. For the summation we 
will do some algebraic transformations to fit the summation formula from Appendix B: 



R {P ,P + \crf\Q + ,Q ) NS = i 



m + -n + ^_-^yi + +m + 



+ X 



Q' + cQ+PLcV + 



x n n ^(«) n °w n -o^Po^KiQ+^fyvs 

pgp; geQ^ peP| geQ' + 

= ^ + -n +( _ l)S++ m +Arp0)Qo J- Yl <Vl,M>< 

Q'+cQ+PtcV+ 



x n m n n 



-a{p) 



tan 2 n ( 

P ev+ q e<2+ pG p; MUV q eQ' + w 



tan — x 
2 



ILe/' (cos x + 1)°" UyGi' ( cos 2/ - !) r iLe/'^el' ( cos x ~ cos 2/) 
where /' = Q' + U Pi, /' = Pi U Qi, 



/i = - (L - m - n ) 



fi = -(L-n - m ) 



The following relations will be useful in what follows: 

m — n 



t + fi = a + fi , |/|=/i 



n — m 



2 ' 1 1 r 2 
where definitions of a and r are given by (|58p . Using the relations 



tan 



P 



/ g(p) 2 - e 2 , a(p) 

4 - 4p) 2 ' Pip) 



tan 



e + e(p) /3(g) 



£tt - e(p) 



a(g) 



tan 



and 



cos x + 1 



4-e(s) ! 



cos x — 1 



e(x) 5 



2h ' 2h 

in the expression for the matrix element we get 



cos x — cos y 



_ £0 ~ e(g) 
Ett + e(<?) 

e(j/) 2 -g(x) 2 
2h 



R (P ,P + \af\Q + ,Q ) NS = i m+ - n+ (-l) n++m++fl N Po>Qo J] /3(p) J] a(g) 

peP+ <?eQ+ 

- e v,,* n f^i n 



p+,Qh 



( 2h )<T|/'|+r|7'|+|/'||J'| 



(73) 



(74) 



(75) 



(76) 
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In the terms of the notations 

7 P = ~e(p), if P G P+> 7g = if 9 € Q+, 

we have 

r{Pq,P+Wi\Q+,Qo)ns = 



i 



m + -n +{ _ ir++ m ++fiNpoQ I ) ] [ ^ - [ a(g)>< 

x V <5 , TT (fZE -^yTico ~ lyf^ 1 

/'c^e+ ^ yei> ^ + 7 ^ )1 " <T (e ° + R*&',ver> hi ~ T# ' 

Now we make summation over I' by means of the following formula from Appendix 

s TT (£tt ~ lyf (gO ~ 7y) 1 " r 1 = 

/'cnuQ + ,i ^ + ho + IvY ' U x er,yeP hi ~ ft) ~ 

2min(/x,/i) /_ jW/x+l)/2 (_]Wl-r) 

rLeP+u2+ + 7*) 1- " (eo + 7*) r IL <2/ eP + uQ + (7* + 7») ' 
where min(/x, /2) = (/i + /2 — |/i — fl\)/2 = + \I\ — \a — t|)/2. As a result we obtain 

R (P ,P + \af\Q + ,Qo)NS = i m+ - n+ (-ir ++fh++p Np , Qo T1 1 J 

Ilx6-P + UQ+ (. e 7T 7xJ 

[ yS(p) [ a(g) 1 ^ , — 

per+ qeQ+ Wx&v+vjQ+hTT + lx) {zq+Ix) Y[ x <yeV + uQ + wa; + 7y) 

Using the relations ([70]) . (f73|) . (|71j) and 7 a = e a for a S /, -yp = —ep for (3 £ I we get 
*<P ,P + ki|Q +) Qo>iVS = i m+ - n+ (^r +Tp+ ^N P(hQo M.x 



{en + e a ){e a + e ) yj j (e n - ep)(ep - e ) 



2e a - V 2efl 

aei v ^g/ V p 

1 



IlaG/ hn + e«) (^tt - £aY (eo + e Q ) T EI/3e7 ( £ 7r - ep) (^ + £pY (~£p + £o) T 

1 

Ua<a'el ( £ " + £ a') Il^'g/ (~ e /3 ~ Dae/ ll/JeJ ( £ « ~ £ /3) ' 

where 

M_ — ^_j\m++M+M(M+l)/2+a-(l-r) _ 
Let us substitute here the coefficient Np 0t Q which follows from (fHTj) of Appendix A: 

/ o \ (n+m)/2 



J m -n o ,0' ^ \r ) 

e ip(m-n)/2 e -ig(m-n-2)/2^^ x 

peP () <?eQo 

(-i) n+ n - !) r ii( cosa + i ) ct n n( c ° sa - ^ . 

06/ /36-T a£/ a£l p<=I 
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where 

A P,Q 



n p < P >ep(e ip -e ip/ )Il q<q >eQ(e iq -e^ 



n ?eQ n pe p(^-^) 

pf — ^_^\n(n— 1)/2/ Wm— n)(m— n— 2)/8/_ mo-no)(n-"i)/4 

Then using ((76j) and 2H< T ~ T I/ 2 = 2" 1 / 4 (e + £7r) (1 ~ 2<T){1 ~ 2r)/4 , which follows from a,r G {0,1} 
and eo = 1 — h, e n = 1 + h, we get 

/ n \ (m+n)/2 

R (P ,P + \af\Q + ,Q ) NS = 5„t°n d o 2) (J) (-l)"^-!)/^.!)^ )/2 h (n-m)V4 x 



X Y[ e ip( - m - n)/2 Yl e -iQ(^-n-2)/2 Ap ^ x 

peP ggQ 



xO-^frn + f \(l-2g)(l-ar)/4 TT / + g Q )( £ « + £ o) TT / i 

n / 3gj(£0+g/3) T (^+g/3) 1 ~ f7 IlagJ^g/ ( g » + £ /?) 



Flag/ ( £ - + z*) 1 ° (eo + e a ) T n Q < Q 'g/ (e Q + £<*') n^'eJ + e /3') ' 

where we also extended the products over Po an d Qo of the exponents of momenta to P and 
Q, respectively, since it means addition of pairs of opposite momenta. Regrouping factors with 
respect to the sets U = Q + P = Q+ + Q- + Q + P+ + P- + P and U = P + Q we obtain 

(r> \ (m+n)/2 
|J (_l)^-l)/2 + (n- m )/2 h (n- m) V4^ QX 

x JJ e *(m-n)/2 JJ e -Mm-n-m (^V r "^'^ n + U/8 " ( 79 ) 

P ei , e Q n Q , a 'gc/(^ + ^) 1/8 iw e &fo + £/3') 1/8 

To prove this formula we have to compare the exponents of (e a +e/3) for all a and f3 in the former 
and latter expressions and to use definitions (i58l) of a and r. Now we rewrite the products over 
the sets U and U as products over the sets P, Q, P, Q. Finally we exclude the products over 
the sets P, Q supplementing them to be the products over all the sets of momenta in the first 
Brillouin zone Q = Q U Q for iVS'-sector and V = P U P for i?-sector. It gives the factorized 
formula for the matrix element of spin operator [4]: 

r{P\ox \Q)ns = r(Pi, P2, ■ ■ ■ ,PnWf\qi, <?2, ■ ■ ■ , 1m) NS = 

= cj r(-+ m )/ 2 (-i) n(n ~ i)/2 n e ~ ip/2 n eiq/2x 

peP <jgQ 
( 2 \ ( m +™)/ 2 . , 2 u l _ e t?(«)/2 e -ij(p)/2 



geQ V -1 pgp V 2e P 



X 



X 



n mh) n S£ n n . <-> 

q < q 'eQ \ £q + £q J P <p'eP \ £ r> + £ P' J geQ peP \ 2sm ~/ 



where e a is given by (|63|) and 

£ = ft _ f, 2 )* ^ = llggQ ]lpgp( £ g + £ p) 2 _ Uq'GQ ( £ « + £ q') 
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Matrix elements with even (resp. odd) n and m correspond to the Hamiltonian (pQ) for periodic 
boundary condition (resp. ([5} for antiperiodic boundary condition). 

Let us comment on the modification of formula (j80|) in the case of the spin <r£ matrix 
elements. From (I68P and similar formula for R-sector and also from (jlip we obtain 

r(P\o%\Q)ns = fl(P|r fe -Vfr- fe+1 |Q)jvs = e^ k ~ 1 ^Q cl ~^pP) R (P\af\Q) N s. 

6 Conclusions 

In the paper [TJ] using Baxter's extension of the Onsager algebra the factorized expressions for 
the spin operator matrix elements between the eigenstates of Hamiltonians of the finite length 
superintegrable ^TV-symmetric chiral Potts quantum chain (SCPS) were found up to unknown 
scalar factors for any pair of the Onsager sectors. In this paper we have derived the exact 
expression (|60|) for these factors for quantum Ising chain in a transverse field (N = 2 SCPC- 
model). This derivation uses some modification of the standard fermion technique [19J. On 
the first stage we diagonalized the Hamiltonian Hq ([5]) by means of the dual Jordan- Wigner 
fermion operators and calculated the spin matrix elements between eigenstates of Hq, then after 
the Bogoliubov transformation we obtained the factorized formula for the spin matrix elements 
between iegenvectors of H using the summation formula (Appendix B). It is natural to expect 
that the fermion technique used in this paper can be applied to calculation of the spin matrix 
elements in more general free fermion models, for example, in the N = 2 Baxter-Bazhanov- 
Stroganov model |20|, [5| . We will address this problem in our forthcoming paper. 

Acknowledgements. We thank G. von Gehlen, O. Lisovyy, S. Pakuliak for useful discussions. 
This work was supported by the Program of Fundamental Research of the Physics and Astron- 
omy Division of the NAS of Ukraine, the Ukrainian FRSF grants $28.2/083 and $29.1/028, by 
French-Ukrainian program Dnipro M17-2009 and the joint project PICS of CNRS and NAS of 
Ukraine. N.I. is thankful to Max-Planck-Institut fur Mathematik for kind hospitality. 



Appendix A 



In this Appendix we prove the following trigonometric identity for the momenta of quantum 
Ising chain. We will use the definitions from Sect. 4.2 for different sets of momenta. Supposing 
that m = \Q\, n = \P\, mo = \Qq\ and no = |Po| have equal parity, we claim that 

, n p<P >ep(e ip -e ip 'm g<q >e Q ^ g -e iq ') 
Ap,Q := 



n geQ n pe p(^-e ip 

\ -(n+m)/2 
_ I „-(no-mo)/- 
J m -n ,0 I £ 



/ o \ ~(n+m)/2 

j( mod2 )/_ | ^(no-m )/2 i (m-n)/2 2 (n-m) 2 /4 ^p{n-m)/2 e iq(m-n-2)/2 x 

pePo q&Qo 



(-l) n+ ELae I tan f 

x N - Np °> Qo rupees a + ir rw^ p - i) r iw(^ « - j) • (81) 

with I = Q + U P+, I = P+ U Q + , 

^(mod2)^ _ 1 ^_^ no - mo ^ N = (— l^ n ( n - 1 )/ 2 +( m - n )( m - n - 2 )/ 8 +( L - m -n )(n-m)/4 

and the coefficient Np 0j q not depending on P + and Q+ (or equivalently not depending on I). 
In the main text we use this formula for different choices of P + and Q + , namely, as labels of 
eigenvectors of Hq or as labels of eigenvectors of H with Pq and Qq being fixed. 
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In order to prove (|8ip it is sufficient to show that addition of a pair of opposite momenta to 
P or to Q in the left-hand side of (|81|) does not change the coefficient Np 0t Q in the right-hand 
side of (|8ip . Let us add, for example, a pair of momenta {— p, p} to P and calculate the ratio 
of left-hand sides of (I8TT) : 



A Pu{ - P , ph Q _ U P ' €Po B(p,p') rW + CW) ( -i P _j P) 

MQ U q 'eQ B(p,q')'U q 'e Q+ C(p,q')' { h 1 j 

where 

B(a, P) = {e ia - e ifi ){e~ ia - e ifi ) , C(a, /?) = B(a, p)B{a, -0) . 

Note that in (|8ip we used the following ordering of the momenta {q\, q2, ■ ■ ■ , q m }- from the 
beginning the pairs (—q, q), for q € Q+, are going on and then q £ Qo in & order fixed in all the 
paper. Analogously we define the ordering of the momenta in the set P = {p±, P2, ■ ■ ■ , p n }- It 
order to present the factors of (|82p in convenient form we introduce the following functions: 



For the products of the functions B we use the first two of following identities: 



(/7_(p) Po ' "^ ' "V+(p) PO '' r ~ PO '' r ~ 1 



n B(P, 9') = 2(-l)™° e-*™° J] eV ^-rf • ( 84 ) 



g'eQo g'eQo 



J] </) = L(-ir e-™ J] ^ ^ n ^ y 7 g) . , (85) 



<?'£Qo <?'e<2o WeS+ 



B(q,p') = 2(-l) no e-^ no TT = . (86) 



p'ePo p'6P 



The functions Pop, Po$, Pq,tt, Po,n are defined by ([59]) . The functions Qo,o, Qo,o, Qo,n, Qo,n 
are defined similarly. Note that always we have Qo,o = Qo,o = but we included them for 
uniformity. The identities (|83p - (|86p follow from the elementary trigonometric identities: 



L-l 



n {is -l, no 



^+1/2-1 = 2 a , = e 2"/l. 



fe = l fc = l 

We also have 

C(a,/3) = e~ 2 * a x(a,(3) 2 , e~*> - e* = e-*<p + (p)tp-(p) 
Since tjiq and no has the same parity, for the ratio (|82p we have 



^PU{~P,P},Q = £ ip(m-n-l) TT e V TT 1-(p) I+Cp) 2 " II/?ejX(P.ff) 

where we used (|58p and the definition of the sets I and J. Taking into account the identities 
e~ ia x(a,P) = 2(cosa - cos/3), e~"V+(a) = 2(cosa + 1), 

(/?) = 2(cos/3 - 1), = -itaa(a/2) 
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and the relation \I\ — \I\ + a — t = n — m following from (|74p we get 



A p\j{-p,p},Q _ L ^n-m+l TT ip' TT 

Apn 2 11 11 

p'GPo q'eQo 



+ p (cosp+iy U^i(cosp-cos^) 

tan - - - ^ . (87 

2 (cosp-l) T n^g(cosp-cosa) 



The corresponding ratio of right-hand sides of (|8ip gives the same result (|87p . For the ratio 
Ap,Qu{-<2,g}/Ap,Q the calculation goes in the same way with the use (|85|) and (|86|) for the 
products of 5. It proves the identity (|8ip . 



Appendix B 

The following summation formula over the subsets / E 1Z, 1Z = V+ U 

>T^^ TT (gyr - lyY (gQ ~ 7y) 1 " T 1 = 

f^ n (e„ + t,) 1 - ( eo + 7s/ ) r ' n.gr.vef (t, 2 - ll) ~ 

2^0*)£)(— 1)m(m+1)/2/_-Q<t(1-t) 

ILeW ( £ t + 7x) 1 ~ <T (go + 7x) r U x <yen (lx + 7y) ' 

where min(/i,/i) = (// + /2 — |/i — /2|)/2, is valid. In fact the formula contains four subcases 
depending on the choice of a,T € {0, 1}. In this formula /i is defined from the relations /i + r = 
/2 + a, [i + p, = \1Z\ = r. It is easy to prove this formula reducing it to the summation formulas 
from Appendix A of [12] for arbitrary variables u, v, z a , a = 1, . . . , r and even and odd r, 
respectively: 

y W (- 1) rlr-2),8 {u _ v y,2 

y *|/|,(r + l)/2 (-l)(H-l)(r-l)/8 (M + t;) (r-l)/2 

We may consider four cases of a, r G {0, 1} one by one. 

For example let us prove (|88p for the case <r = 1, r = 0. In this case jjL = (r + l)/2, 
/2 = (r - l)/2. The left-hand side of {88} is 

ICR ilx€l,y€l \'y IxJ y£ Ti 
^ <Vl,(r+l)/2 (-1) MA = 2( r - 1 )/ 2 (-l)^-l)/2+/x-l 

n^e/ (7y - gyr) (7j/ - go) U x ei,yel bl ~ ly) Ux< y en (7* + 7») 

where we used (f90|) for u = — eq, v = —e n , {z a } = {7a} and u + v = —2. Thus we proved 
for the case a = 1, r = 0. The other three cases of a and r can be considered similarly. 
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